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ABSTRACT 

We measure the strength, frequency, and timescale of tidally triggered star formation 
at redshift z — 0.08 — 0.38 in a spectroscopically complete sample of galaxy pairs drawn 
from the magnitude-limited redshift survey of 9,825 Smithsonian Hectospec Lensing Survey 
(SHELS) galaxies with R < 20.3. To examine the evidence for tidal triggering, we identify a 
volume-limited sample of major (|AMr| < 1.75, corresponding to mass ratio > 1/5) pair galaxies 
with Mr < —20.8 in the redshift range z — 0.08 — 0.31. The size and completeness of the 
spectroscopic survey allows us to focus on regions of low local density. The spectrophotometric 
calibration enables the use of the the 4000 A break (D„4000), the Ha specific star formation 
rate (SSFRy Q ), and population models to characterize the galaxies. We show that D„4000 is a 
useful population classification tool; it closely tracks the identification of emission line galaxies. 
The sample of major pair galaxies in regions of low local density with low D„4000 demonstrates 
the expected anti-correlation between pair-wise projected separation and a set of star formation 
indicators explored in previous studies. We measure the frequency of triggered star formation 
by comparing the SSFRh^ in the volume-limited sample in regions of low local density: 32 ± 7% 
of the major pair galaxies have SSFR# Q at least double the median rate of the unpaired field 
galaxies. Comparison of stellar population models for pair and for unpaired field galaxies implies 
a timescale for triggered star formation of ~ 300 — 400 Myr. 



Subject headings: galaxies: interactions - galaxies: 
1. Introduction 

Hierarchical galaxy formation models predict 
frequent interac t ions in a galaxy's h istory (e.g. 
Cole et all I2OO0I: IWechsler et al.l (20021) . Some of 
these interactions may trigger star formation. Ob- 
servational limits on the frequency and duration 



stellar content - galaxies: active 

of triggered star formation are an important con- 
straint on the role of these interactions in the evo- 
lution of stellar populations. 

Evidence that galaxy interactions at low red- 
shift trigger star formation was first detected 
in photometric observations by lLarson fc Tinslev 
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(|l978l ), who demonstrated that apparently inter- 
acti ng galaxies in the Atlas of Peculiar Galax- 
ies |a3 Il966l) have a broader range of colors 
and tidal features than typical galaxies. Spectro- 
scopic indicators confirm that interacting galax- 
i es tend to have enhanc e d star formation rates 
Kennicutt fe Keell Il984t iKennicutt et all Il987t 



Keel 



1997 



2003 



1993tlLiu k Kennicuttlll995a| : lDonzelli k Pastorizal 



Barton et al 



Nikolic et al. 



2008) . Infrared observations (IKennicutt et al 



2000, 



■2003t Lambas et al 
2004 iKauffmann et al.l 12004 



Woods et alf2 006; Wo ods fc Gellerll2007l;lElhson et al 



1987 : 1 Jones fe SteiiJll98ilSekiguchi fe Wolstencroft 



1992 : iKeell Il993t iNikolic et all 12004 lOeller et al 
20061: ISmith et al.l l2007) and radio observations 
(jHummell Il98ll) yield similar results. Observa- 
tions that galaxy pairs with the smallest project 
separation have the highest star formation rates 
provide direct evidence for triggered star forma- 
tion (e.g Barton et al. 2000l) . 

Recent studies show that galaxy interactions at 
intermediate redshift trigger star formation. Mea- 
surements of infr ared luminosity dLin et alJl2007h 
and EW([0 II]) (|de Ravel et alJl2009h indicate an 
anti-correlation between projected separation and 
star formation activity for pair galaxies at inter- 
mediate redshift (z < 1), similar to the trend ob- 
served at low redshift. Analysis of the morphol- 
ogy of intermediate redshift (z = 0.1 — 0.6) pair 
galaxies likewise shows an increase in the fraction 
of asymmetric g alaxies in pairs, an indication of 
tidal interaction ( Patton et al. 20051) . 

Although observations demonstrate that galaxy 
interactions can lead to enhanced star formation 
activity, the results of studies of the frequency and 
intensity of star formation a ctivity appear to vary 
with the sample selection. iLi et al.l ( 2008a ) find 
that 30% of the SDSS pair galaxies with high 
specific star formation rates have a companion 
within 100 kpc, and that low to average specific 
star formation rate galaxies rarely have a com- 
panion. In a se parate study of SDS S pair galax- 
ies at z < 0.16, Ellison et al. I (|2008h measure up 
to 40% enhancement in the median star forma- 
tion rate of close pair galaxies with stellar mass 
ratios < 1/10 relative to their comparison sam- 
ple of non-pair systems. iBergvall et al. I (l2003h find 
that global U BV colors do not show significant en- 
hancement in star formation in their small sample 
of interacting galaxies compared to isolated sys- 



tems. However, they do find that star formation 
rates at the very centers of the interacting sys- 
tems are increased by a factor of ~ 2 — 3 over the 
non-interacting systems. Spitzer infrared obser- 
vations of 35 tidally disturbed Arp galaxies show 
that their 24 /iin emission is more centrally con- 
centrated than in normal spiral galaxies, suggest- 
ing a build up of central gas which could fuel cen- 
tral star formation; their infrared colors suggest an 
increase in the mass normaliz ed SFR by a facto r 
of two over the normal spirals (jSmith et al.ll2007l) . 

Numerical simulations account for the observed 
range of strength, frequency, and duration of trig- 
gered star formation. In the larg e suite of sim- 
ulations bv iDi Matteo et al.l (|2008h . galaxy inter- 
actions and mergers produce only moderate star 
formation enhancement relative to non-interacting 
galaxies at low redshift. Strong starbursts are rare 
and short lived, typically lasting a few hundred 
Myr. The duration of the starburst declines as 
the enhancement in star formation rate increases. 

Global galaxy properties affect their suscepti- 
bility to tidally triggered star formation. Systems 
with sufficient gas can exhibit t riggered star for- 
mation in major interactions (e.g. Mihos fe Hernquisd 



1996tlTissera et al.l2002llCox et alj2006HDi Matteo et al 



120071 ). Interactions between gas-poor galaxies 
("dry mergers") produce little or no star forma- 
tion, although they contribute su bstantially to the 
build-up of massive galaxies (e.g. Tran et alj|2005t 



Van Dokkum 



issive g 
20051 : K 



Cattaneo et al 



20081) . 



Inter- 
nal structure strongly affects the susceptibility to 
gaseous inflows; late- type galaxies without strong 
bulge components are more likely to have bar 
instabilities that drive the gas eous inflows power- 
ing the central star formatio n ()Mihos k Hernquist 
19961: IDi Matteo et al.ll2007h . 



The relative mass of the galaxy compared to 
its neighbor and the intrinsic mass also influence 
the degree of triggered star formation. Galaxy 
pairs of similar mass (luminosity) show strong 
enhancement to their star formation rates in 
both members of the pair in both observa tional 
(IWoods fc Gelled l2007t lEllison et al.1 120081 ) and 
theoretial ( Cox et all 20081 ) studies. Observations 
show that although satellite galaxies occasionally 
experience triggered star formation, the brighter 
of a high mass (luminosity) ratio pair does not un- 
dergo significant tidal triggering ([Woods fc Geller 
20071 ) ; numerical simulations suggest a similar pic- 
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ture ( Cox et al. 20081) . The intrinsic mass of the 
galaxy also plays a role: low mass (luminosity) 
galaxies in major mergers are more strongly af- 
fected by the interaction than are high mass galax- 
ies i n major mergers in ob s ervations of pair ga lax- 
ies (IWoods fc Gellerll2007l [Ellison et al.ll2008h . 

Local environment can affect star formation 
activity in a galaxy apparently independent of 
an interaction. Galaxies in clusters and large 
groups tend to have redder color and less cur- 
rent star formation than is olated system s (e.g . 
Hubble fe Humasonl Il93l ICooper et all 120071 
Gerke et al.ll2007l ). iBarton et all ([20071 ) show that 
pair galaxies are more common in regions with 
greater local density than are unpaired galaxies. 
Including pair galaxies from high density regions 
suppresses the signal of triggered star formation 
measured in comparison to unpaired galaxies. 

We take advantage of the size and complete- 
ness of the Smithsonian Hectospec Lensing Sur- 
vey (SHELS. iGeller et al.|[2005l l2010h galaxy sam- 
ple to identify a set of major interactions at z = 
0.08 — 0.38 for galaxies with young stellar popu- 
lations in regions of low local density, where we 
can measure the cleanest signal of triggered star 
formation. Local density selection at intermediate 
redshift has not been possible before this sample. 

The stability of the Hectospec instrument en- 
able accurate spectrophotometry calibration of 
the spectra. The calibrated spectra provide the 
4000 A break (£>„4000) and the Ha specific star 
formation rate (SSFRhq), which we use to char- 
acterize the star formation. We quantify the fre- 
quency, strength, and timescale of triggered star 
formation in major galaxy interactions using mea- 
surement of specific star formation rates and stel- 
lar population models fitted to the individual spec- 
tra. 

We describe the SHELS data in jJU and the 
sample selection in Section 0] describes the 
classification of galaxies as starforming or AGN 
and the use of the spectroscopic indicator Z?„4000 
to identify galaxy types. We consider local den- 
sity effects in £j5j Section [6] describes various 
measurements of tidally triggered star formation, 
including the frequency, strength, and timescale 
of enhanced star formation activity. Throughout 
this study we assume standard ACDM cosmology, 
where Hp = 71 km s" 1 Mpc -1 , f2 m — 0.3, and 



2. Data 

In this section we describe the Smithsonian 
Hectospec Lensing Survey data set (SHELS). We 
give an overview of the data set in ^2.11 We de- 
scribe the photometric measurements in £12.21 and 
the spectroscopic measurements in §2.3[ including 
the Ha specific star formation rate in §2.3. 1[ and 
aperture effects in i j2.3.2l 

2.1. Smithsonian Hectospec Lensing Sur- 
vey 

The SHELS magnitude-limited redshift survey 
includes 9,825 galaxies with total i?-band magni- 
tude R < 20.3. The galaxies are selected from the 
Deep Lens Survey F2 field R-band images (DLS; 



Wittman et all I2002L 120061). We obtained spec 



troscopy using Hectospec ()Fabricant et al.l [20051 ) 
at the MMT. The spectroscopic sample is 97.7% 
complete to 11=20.3, and the differential complete- 
ness at the magnitude limit is 94.6%. 

We restrict our analysis to the 6,935 galaxies 
with redshift z = 0.080-0.376. The lower redshift 
limit minimizes aperture effects ( H2.3.2p and the 
upper limit allows measurement of the Ha flux 
(» 

The DLS F2 field covers 4 deg 2 on the sky cen- 
tered on R.A. = 9 h 19 m 32 s 4, decl.=+30°00'00" 
(J2000). The Sloan Digital Sk y Survey (SDSS 



Adelman-McCarthv et al. 20061 ) also covers this 



0\ = 0.7 (|Spergel et al.l 120031 ). 



region. 

The DLS images come from the Kitt Peak May- 
all 4-meter telescope with the Mosaic prime-focus 
imager (Muller et al. 1998). The DLS observed 
in the i?-band during nights with seeing < 0'.'9. 
The 1-er limiting surface brightness in i?-band is 
29 mag arcsec -2 . The DLS also observed in V- 
band. 

The 0"9 resolution of the DLS i?-band im- 
ages allows identification of close pairs or merg- 
ing galaxies with minimum separations of 15 — 
20 kpc up to z < 0.3. The median seeing 
for the SDSS photometry is 1" 4 PSF in r-band 
( Adelman-McCarthv et al. 20061 ). Thus some sys- 
tems that cannot be distinguished in the SDSS 
data can be identified in the DLS. A total of eight 
pairs in the DLS are missing from SDSS. 

Objects with stellar light profiles, i.e. AGN 
at high redshift, may be preferentially excluded 
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from our sample because DLS objects with stel- 
lar light profiles are not targeted for spectroscopy. 
However, with the DLS 29 mag arcsec -2 surface 
brightness limit and the 0'.'9 resolution, we should 
be able to detect bright host galaxy bulges to the 
limiting z — 0.376. We discuss AGN detection in 
more detail in ij4.ll 



The Hectospec observations identify a spectro- 
scopic pair of emission line galaxies at R.A.= 
9 h 16 m 58 s 903, decl. = +29°43'10'/597 (J2000). 
The image, which shows a double-nuclei extended 
object, and the spectra are available on our web- 
sitdE We do not include this object in our analysis 
because the selection of spectroscopic pairs is not 
complete or uniform. 

Only 230 galaxies in our sample lack spectra 
(2.3%). Of these, - 24 galaxies (10%) may be 
members of pairs that satisfy our criteria for anal- 
ysis (|J3|). The fraction of galaxies in pairs in our 
sample as a whole is ~ 10% (!}3J). 

2.2. Photometric measurements 

Image processing for the DLS uses the IRAF 

package MSCRED to correct for bias, to fla t field, 

and to obtain basic astrometric calibration. Wittman et al 



photometry because the shape of the filter band 
passes are similar and the difference in corrections 
for the different bands is negligible. We estimate 
the uncertainty in the k+e correction from the 
difference in the correction for adjacent galaxy 
types. The uncertainty ranges from a maximum 
of ~ 0.02 mag at z = 0.08 to 0.03 - 0.13 mag at 
z = 0.38, depending on the galaxy type. 

There are 51 galaxies in our z = 0.080 — 0.376 
sample with a Hectospec redshift but no SDSS 
photometry for galaxy classification and k+e cor- 
rection, including three galaxies in major pairs 
(jpT2]l, In 34 of the 51 cases, we use the DLS 
(V — R) color to estimate the (g — i) color for the 
k+e correction. For the 17 cases where we lack 
both DLS (V - R) and SDSS data (including two 
in major pairs), we assume the galaxy type is Sa. 
The k+c correction for the galaxy type Sa is in the 
middle of the range for the different galaxy types; 
choosing a different k+e correction has no signifi- 
cant effect on the results. At the redshift z = 0.3, 
the apparent R — 20.3 limit corresponds to a k+e 
corrected absolute magnitude Mr = —21.4 for a 
typical late- type galaxy, and to Mr — —21.8 for a 
ypical early-type galaxy. 



(|2006h construct a stacked DLS image of each 
subfield in the i?-band, correcting for cosmic rays 
and saturated objects, fixing the astrometric cal- 
ibration, and correcting the image shape and 
photometry for optical distortion. The uncer- 
tainty in the i?-band magnitude is ~ 0.05 mag. 
We derive the galaxy c atalog using SExtractor 
(jBertin fe Arnoutslll996f ) on the final stacked im- 
ages. We remove objects within the radius of the 
diffraction pattern around bright stars. The total 
area excluded around 778 stars is 0.215 deg 2 (5.4% 
of the survey). 

The absolute magnitude Mr calculation in- 
cludes both k- and evolutionary (e) corre ctions. 
We us e the k+e correction determined by Annis 



2.3. Spectroscopic Measurements 

The Hectospec spectrograph 270 line mm -1 
grating yields ~ 6 A pixel -1 dis persion over the 



wave length range 3700-9300 A (Fa bricant et al 



20081 ). The 300 optical fibers with l'/5 diame- 
ter are placed radially within a 1° diameter field. 
About 30 of the fibers measure sky background 
during each pointing. 

We reduce Hectospec data with the standard 
Harvard-Smithsonian Center fo r Astrophysics 
Hectospec reduction p ipelinjl (jFabricant et al 



2005: iMink et all 120071 ). We meas ure the galaxy 



redshift using the program RVSA0 ( Kurtz fe Mink 
ll998tlFabricant et al.ll2005T ). The total RMS inter- 
(l200ll) using the Pegase code (|Le Borgne fc Rocca-Volmerag# error in the redshift is 34 km/s for emission line 
20021 ). The k+e correction requires classification galaxies and 65 km/s for absorption line galaxies, 

as one of nine galaxy types (bright cluster galaxy, based on repeat observations of 812 emission lines 

elliptical, SO, Sa, Sb, Sbc, Sc, Sd, and irregular). and 542 absorption line galaxies. We further pro- 

We classify the galaxies according to the SDSS cess spectra to correct for atmospheric extinction 

(g - i) color and galaxy redshift. We apply the and to remove the narrow absorption lines from 

SDSS r-band k+e correction to the DLS i?-band H 20 and 2 at wavelengths longer than - 6000 A. 



1 http://tdc-www.cfa.harvard.edu/instruments/hectospec/progs/tsf/ 2 http://tdc- www. harvard.edu/instruments/hectospec/reduce.html 
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We also correct for relative throughput a s a func- 
tion of wavelength ( Fabricant et al. 20081 ) . 



We compute the emission line flux by integrat- 
ing the counts in a fixed-width band in the rest- 
frame centered on the central wavelength of the 
line; we subtract the continuum level. The spec- 
tral regions where we measure Ha and H/3 emis- 
sion are 6562. 8±8 A and 4861. 3±8 A, respectively. 
We compute the continuum from the average of 
regions on either side of the band after rejecting 
noisy data. The RMS scatter in the line flux is 
23% for Ha, and internal systematic error in Ha 
is ~ 18%, based on duplicate observations of 592 
galaxies. 

The calibrated i?-band photometry from the 
DLS in a l'/5 diameter aperture serves as reference 
for conversion of counts to flux in the l'/5 diam- 
eter fiber. This calibration process removes vari- 
ation caused by seeing, clouds, telescope tracking 
and guiding, and astrometry and alignment errors 
for light losses that a re wavelength independent 
(jFabricant et al.ll2008t ). 

Comparison of Hectospec spectra with a sam- 
ple of overlapping SDSS spectra shows that the 
median ratio of the spectra is constant to ~ 5% 
over the range 3850 — 8000 A, and the offsets in 
the ratios of the median Ha flux and median rest- 
frame Ha equivalent w idth (EW(Ha)) are ~ 3% 
(jFabricant et al.ll2008t ). 

We also compute the spectroscopic indicator 
D„4000, a measure of the stellar population age, 
which is defined as the ratio of flu x in the band 
4000- 4100 A to flux in 3850-3950 A (|Balogh et al 
19991 ). The RMS scatter in Z?„4000 is 0.086 



The internal systematic errors in D n 4000 are very 
small, 4.5%, based on repeat measurements of 
1468 galaxies. The comparison of our measured 
D„4000 with an overlapping samp le of SDSS spec 



tra shows a median ratio of 1.00 ( Fabricant et al 



2008). 



2.3.1. Ha specific star formation rate 

We compute the Ha specific star formation rate 
(SSFR# a ) from the Ha flux (fija) normalized by 
the i?-band 1'.'5 aperture magnitude. The Ha flux 
must first be corrected for stellar absorption and 
for reddening. 

To determine the stellar absorption at Ha 
and H/3, we fit the individual spectrum using 



tion method (see also IWestra et al 
method fits the calibrated spectra 



the Tremonti et al. (2004) continuum subtrac- 
tion This 
with stel- 
lar population models to estimate the ages of 
the stellar po pulations within the galaxy. Like 



Tremonti et all we assume that the star forma- 
tion history of a galaxy can be represented as the 
sum of discrete bursts of star formation, and we 
fit for the burst ages. 

Figure [1] illustrates the fit of the contributions 
from the different starburst populations in an ex- 
ample galaxy s pectrum. The template spectra 



that we use are iBruzual &: Charlotl (|2003l ) stellar 
population models. The models include 10 dif- 
ferent ages of bursts (0.005, 0.025, 0.1, 0.3, 0.6, 
0.9, 1.4, 2.5, 5 and 10 Gyr) at Solar metallicity. 
We also add r eddening as an additional param- 
eter using the ICharlot fc Falll (|2000j) attenuation 
model. Finally, we match the models in redshift, 
pixel size, and spectral resolution to each galaxy 
spectrum. Th e core of the code adapted from 



Tremonti et al 



uses the ID L mpfitfun rout ine 
from the Markwardt library ( Markw ardt 200jJ to 
find the best fitting combination of models. 



1 40.9830938 + 30.691 4996 




5000 5500 
Wavelength (A) 



Fig. 1. — An example rest- frame spectrum illus- 
trating the contribution of different stellar popula- 
tions to the continuum galaxy spectrum. The solid 
black line is the observed galaxy spectrum. The 
solid red line shows the modeled galaxy continuum 
from the sum of the discrete starburst populations, 
shown as dashed and dash-dot lines (blue). 
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The absorption-corrected Ha flux is: 
(EW{Ha) + EW Ah , 



/ Cl Ha f Ha y EW(Ha) 



(1) 



and similarly for / c (H/3). The ratio of / c // is 
equivalent to EW C /EW. All EW and line fluxes of 
Ha or H/3 in this paper are absorption corrected. 

We correct the Ha flux for reddening accord - 
ing to the standard method ( Calzetti et al. 2000). 
We compute the wavelength dependent extinc- 
tion, A\, where the factors k(Ha) — 3.326 and 
k(H/3) = 4.598 give the differential extinction be- 
tween the wavelengths of Ha and H/3 for the case 



of a s tarburst galaxy (k(V) = 4.05. ICalzetti et aT 



2000). The parameter i? Q( g relates the observed to 
intrinsic Balmer decrement (2.87; Calzettill200ll ). 

We correct for reddening using average E(B — 
V) values for bins of absolute luminosity and red- 
shift. This procedure mitigates the effects of noisy 
individual H/3 measurements, which contribute 
significant noise to the E(B — V) values. Fig- 
ure[2]shows the EW(Ha) versus Balmer decrement 
for individual galaxies. A population of galaxies 
with EW(Hq) < 20 A and Balmer decrement 
> 10 result from noisy H/3 measurements. The 
total internal error in EW(Ha) is 18% and the 
error in the Balmer decrement is 25%. We deter- 
mine the Balmer decrement for each luminosity 
and redshift bin from the median value in the ap- 
propriate bin derived from the 484 galaxies with 
S/N> 3 in Ha, S/N> 2 in H/3, EW(Ha) > 3 A, 
EW(H^) > 2 A, D„4000 < 1.44, and classified as 
starforming fi )4.1|) . Using the median E(B — V) 
value eliminates huge and unphysical reddening 
correction factors for individual galaxies. 

Once we correct the Ha flux for stellar absorp- 
tion and for reddening, we calculate the Ha lumi- 
nosity (Lffa) using the luminosity distance. We 
convert from Ljj a t o star formation r ate using the 
conversion factor of Kennicutt ( 19981) : 



SFR Ha (M Q yr- 1 ) 



■>Ho 



1.27 x 10 41 erg s" 1 



(2) 



We normalize the Ha star formation rate by the 
galaxy i?-band aperture magnitude to compute 
the Ha specific star formation rate (SSFR# a ). We 
use the l'/5 diameter i?-band aperture magnitude 
from the DLS to match the fiber diameter where 
the Ha flux is measured. Thus, no correction for 




50 100 

EW(Ha) (Ang) 



200 



Fig. 2. — EW(Hq) versus Balmer decrement for 
all emission line galaxies. The galaxies with 
EW(Ha) < 20 A and Balmer decrement > 10 
result from noisy H/3 measurements. We use the 
median Balmer decrement in bins of absolute lumi- 
nosity and redshift. We set all Balmer decrement 
< 2.86 to that value. Measurement uncertainty in 
EW(Ha) is 18% and the uncertainty in the Balmer 
decrement is 25%. 



aperture covering fraction is needed. The normal- 
ization factor assumes I Mq = lL^, in the i?-band 



and M R{Q) = 4.42 (|Binnev fc Merrifieldl[l998h 



The uncertainty in the SSFR# Q follows from 
the 23% RMS error in the Ha flux and the 
0.05 mag uncertainty in the normalization. The 
absorption and reddening corrections contribute 
additional uncertainty of ~ 20%. We add these 
contributions to the error in quadrature to esti- 
mate a total uncertainty of ~ 30% in estimates of 
SSFRh 

2.3.2. Aperture effects 

The l'/5 Hectospec aperture implies that the 
aperture covering fraction varies as a function 



of galaxy luminosity and redshift. Kcwlcy et al 



( 20051 ) suggest that aperture covering fractions 
greater than 20% yield star formation rates 
representative of the galaxy as a whole. The 
aperture size required to meet this standard 
varies with both redshift and absolute luminos- 
i ty because brighter galaxies are generally larger 
( Kewlev et al. 1 12005 ). In the volume- limited sam- 
ple (defined in i j3.2p . the galaxy radius included 
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in the fiber varies from r = 1.1 kpc at the lower 
redshift limit {z = 0.080) to r = 3.4 kpc at the 
upper redshift limit of the volume-limited sample 
(.2 = 0.310). 

We compute the covering fraction from the dif- 
ference between the 1'.'5 aperture magnitude and 
the total galaxy magnitude. Most galaxies in 
the sample have a covering fraction > 20%; Fig- 
ure [3] shows the fraction of galaxies with cov- 
ering fraction > 20% as a function of redshift 
for the galaxies (both paired and unpaired) at 
z = 0.080 - 0.376. 




,8 0.0 



Fig. 3. — The fraction of galaxies where the l'/5 
aperture includes > 20% of the galaxy light. 



Fabricant et al. (l2008h conduct a detailed in- 



vestigation of aperture effects in their sample of 
overlapping Hectospec and SDSS spectra, a sub- 
set of our full spectroscopic survey. SDSS spec- 
tra are acquire d through a 3" aperture. Thus 
Fabricant et al.l scale the Hectospec line fluxes by 
the 3" aperture DLS i?-band magnitude to make 
direct comparisons with the SDSS line fluxes. 
They conclude that the median [O II] 3727 A and 
Ha line fluxes from Hectospec agree to within 
3 — 4% of the SDSS values, and that the scat- 
ter in the line fluxes is dominated by systematic 
uncertainties. The 3 — 4% difference between Hec- 
tospec and SDSS line flux measurements is small 
compared to the 23% RMS scattter in the Ha line 
flux from Hectospec. 

The agreement between the Hectospec and 
SDSS data implies that the Ha star formation 
rate within apertures of 1.5 and 3" scale with the 
i?-band galaxy luminosity. The relative agree- 



ment between the scaled Hectospec and SDSS Ha 
line flux is indepe ndent of the absolute magnitude 
(see Figure 13 in lFabricant et al.ll2008l ). apparent 
magnitude, redshift, and 25% light radius of the 
galaxies. Westra et al. (2009, in prep.) similarly 
find that the Ha luminosity determined from the 
1'.'5 fiber and i?-band magnitude is consistent with 
the Ha luminosity from the 3 " aperture photome - 
try and i?-band magnitude of lShioya et al.l (|2008l ) . 

Like the Ha luminosity, we assume that D„4000 
observed in th e l'/5 aperture re p resent s the galaxy 
as a whole. iFabricant et al. ( 2008 ) show that 
the Hectospec and SDSS measurement of -D„4000 
agree very well over the scales observed (1.4 to 
2.8 kpc radius at z = 0.1). This agreement again 
suggests that measurement of D„4000 for the stel- 
lar population included in the Hectospec aperture 
is not significantly affected by systematic biases 
compared to the population observed in the larger 
SDSS aperture. 

3. Pair and Field Sample Selection 

We measure the effects of galaxy interactions 
on the star formation activity and other galaxy 
properties at intermediate redshift. We use two 
pair samples derived from the 6,935 galaxies with 
spectroscopy and R < 20.3 in the redshift range 
z = 0.080 - 0.376. 

The first sample, the "full" major pairs sample, 
includes all galaxies in major pairs that meet our 
projected spatial and line-of sight peculiar veloc- 
ity criteria ( 33~Tj) in the range z = 0.080 - 0.376. 
There are 622 galaxies in the full major pairs sam- 
ple. 

The second sample, the "volume-limited" ma- 
jor pairs sample, is a subset of the full major pairs 
sample, and is restricted to galaxies that meet 
both redshift and absolute luminosity selection cri- 
teria ( fl3.2j) . The volume- limited major pairs sam- 
ple allows us to study the strength, frequency, and 
timescale of triggered star formation. The volume- 
limited major pairs sample is well suited to com- 
parison wit h the predictions of the numerical sim- 
ulations of iDi Matteo et al.1 (|2008l ). who measure 
the intensity, frequency, and duration of merger- 
driven star formation in their large suites of nu- 
merical simulations of major interactions. The 
volume-limited major pairs sample contains 339 
galaxies. 
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Our pair samples are the largest spectroscopic 
samples of pairs to date at intermediate redshift 
(z ~ 0.1 — 0.4). In compar ison, the Millen- 
nium Galaxy Cat a log fMGC: iLiske et alj 12003 ; 
iDriver et all 120051: I Allen et al.1 l2006h is a 96% 
complete spectroscopic survey of 10,095 galax- 
ies to B M GC < 20 mag. The MGC is sim- 
ilar in number and completen ess to our sam- 
ple b ut is at lower redshift. De Propris et al.l 
( 2007t ) study pair galaxies in the MGC: they 
identify a volume-limited sample of 3,237 galax- 
ies in the range 0.010 < z < 0.123, including 
112 galaxies in pairs with projected separation 
AD < 28 kpc [20 h^ kpc}. The CNOC2 Redshift 
Survey contains redshifts for ~ 5,000 galaxies at 
0.1 < z < 0.6 with a cumulative completeness 
of 50% and differe ntial completeness of 20% at 
the limit R„ < 21.5 lYee et al.ll200ol: IPatton et aT 



20021 ). IPatton et al.l (|2005l ) study properties of 



dynamically close pairs in the CNOC2 survey. 
Our sample covers the same redshift range but is 
substantially more complete than CNOC2. Ta- 
ble [1] lists the selection criteria and the number of 
galaxies in the full major pair sample and in the 
volume-limited major pair sample. 

3.1. Full major pairs sample 

We select galaxy pairs with projected spatial 
separation AD < 70 kpc, corresponding to the 
limit AD < 50 hTnn kpc commonly appl ied in pre- 
vious studies (e.g. iBarton et ail l2000h . We re- 
quire a line-of-sight peculiar velocity difference 
AV/(1 + z) < 500 km s" 1 . The limit of AD < 
70 kpc includes pair galaxies with the most signifi - 
cantly enhanced EW(Ha) fe.g lBarton et all2 000). 
while minimizing the presen ce of interlope r s. Th e 
AV limit is motivated by IBarton et all (|2000l ). 
who find no significantly enhanced EW(Ha) emis- 
sion in pairs with AV > 500 km s _1 (their Fig. 
2b) in their sample of 786 galaxies from the CfA2 
Redshift Sur vey, and is consistent with previous 
pair studies (|De Propris et al. 2007 ; Woods et al 
20061 ). 



We exclude potential pair galaxies from our 
sample if either member of the pair lies within a 
projected 70 kpc from the survey boundary or the 
edge of the region excluded around a bright star. 
For galaxies in a compact group, we compare each 
galaxy to its nearest neighbor; odd numbers of 
"pair" galaxies can occur in compact groups. 



We compute the total pair fraction for all pairs 
- major and minor - in the redshift interval z = 
0.080 — 0.376. There are 809 galaxies in major or 
minor pairs or compact groups that meet the pro- 
jected spatial and peculiar velocity requirements. 
The total pair fraction is 12% (809/6,935). 

We restrict our sample to major pairs with 
|AA/r| < 1.75, corresponding to a luminosity 
(approximate mass) ratio > 1/5. Simulations of 
interacting galaxies predict that major interac- 
tions are more effective than minor interactions 
at triggering star formation; only the fainter com- 
panion in a minor interac tion occasionally shows 
triggered star formation ( Cox et al.l 120081 ). Ob- 



servations of major and mino r pair galaxies sup- 
port this predicted behavior (Dasvra et al. 20061: 
Woods fc Gelien 120071 : lEllison et al.1 120081 : Westra 
et al. 2009, in preparation.) The full major pair 
sample includes 622 galaxies in major pairs in the 
range z = 0.080 — 0.376. Figure 0] shows the red- 
shift distribution and Figure [5] shows the distri- 
bution of projected separation for the full major 
pairs sample. Table[2]lists the full major pair sam- 
ple galaxies and their properties. 



£ 40 




0.10 0.15 0.20 0.25 0.30 0.35 



Fig. 4. — Redshift distribution for the full major 
pair galaxies (open) and for the volume-limited 
major pair galaxies (filled). 



Visible evidence of tidal interactions in the DLS 
combined BVR-ba,nd images, such as tidal tails, 
marked asymmetries, or extended halos, can be 
seen in ~ 1/4 of the major pair galaxies. Analyz- 
ing the morphological features in detail is beyond 
the scope of this project. The images are available 
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Table 1 

Criteria and size of galaxy pair samples. 



Criteria 


Full sample 


Volume-limited sample 


Major pairs a 


622 


339 


Non-AGN b 


601 


324 


£> n 4000 < 1.44 c 


302 


134 


N c < 4 d 


214 


72 







Note.- The sample size includes galaxies that meet 
listed condition and all previously listed conditions. 

*\AM R \ < 1.75. 

b AGN classification in gUJ 

C D„4000 classification in U7Z[ 

d Local density measure N c described in $5] 



Table 2 

Sample of major galaxy PAiRS a ' fc 



Pair 


RA C 


Dec c 


z d 


m R d 


M R d - e 


AD 




£>„4000 d 


SSFRjj a d 














(kpc) 


(km s _1 ) 




(M yr-yi0 10 Mq) 


001 


9:14:52.622 


30:05:40.858 


0.26472 


19.56 


-21.23 


68 


388 


1.29 


1.36 


001 


9:14:53.920 


30:05:40.377 


0.26309 


19.97 


-20.76 


68 


388 


1.79 


0.10 


002 


9:14:52.849 


30:17:05.567 


0.12805 


20.10 


-18.80 


50 


51 


0.96 


1.26 


002 


9:14:53.528 


30:17:25.933 


0.12813 


19.17 


-19.74 


50 


51 


1.05 


3.99 


003 


9:14:53.583 


30:31:44.906 


0.25643 


20.20 


-20.40 


65 


108 


1.14 


1.18 


003 


9:14:54.188 


30:31:30.369 


0.25688 


19.60 


-21.14 


65 


108 


0.00 


0.04 


004 


9:14:57.683 


29:57:15.864 


0.24721 


19.08 


-21.44 


26 


133 


1.07 


3.83 


004 


9:14:58.198 


29:57:15.988 


0.24780 


19.21 


-21.34 


26 


133 


1.11 


3.79 


005 


9:14:57.762 


29:49:34.947 


0.18074 


18.58 


-21.20 


26 


171 


1.44 


0.36 


005 


9:14:58.434 


29:49:34.583 


0.18014 


17.81 


-22.01 


26 


171 


1.73 


0.00 



a Full table available online: j http: / /tdc- www. cfa.harvard.edu/instruments/hectospec / progs / tsf/| 

b Standard ACDM cosmology: H = 71 km s" 1 Mpc -1 , fi m = 0.3, and Q\ = 0.7. 

C J2000 coordinates. 

d Typical error estimates: 
cz : 34 km s _1 emission line galaxies, 65 km s _1 absorption line galaxies f ij2.3l) . 
m R : 0.05 mag f ^2~2l) . 
£> n 4000 : 4.5%. 

Mfl :~ 0.1 mag, varies with redshift and galaxy type fi j2.2l) . 
SSFR ffQ : 30% ( CI) . 

e k+e corrected. 
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on our website . 

3.2. Volume-limited major pairs sample 

We construct a volume-limited sample of major 
galaxy pairs to study trends in triggered star for- 
mation and other galaxy properties. The volume- 
limited sample contains galaxies in the redshift 
range 0.080 < z < 0.310 and with magnitude 
Mr < —20.8. The limits on Mr and z maximize 
the number of galaxy pairs in the volume-limited 
sample. The lower redshift limit, z = 0.080, min- 
imizes aperture effects fi j2.3.2[) . There are 339 
galaxies in major pairs in the volume-limited sam- 
ple. 

The distribution of projected separations (Fig- 
ure |SJ) is relatively flat from 15 < AD < 70 kpc, as 
expected from the g alaxy correlation function (e.g. 



Estrada et al. 20081 ). At separations < 15 kpc, 



galaxies in the DLS images may not be resolved 
into separate systems. Figures |6] shows the distri- 
bution of absolute magnitudes of galaxies in the 
full major pair sample and in the volume-limited 
major pair sample. 




30 40 
AD (kpc) 



Fig. 5. — Distributions of projected separation for 
galaxies in the full major pair sample (open) and 
for the volume-limited major pair sample (filled). 
The distribution is nearly flat from 10 < AD < 
70 kpc, as expected. At very small AD, not all 
pairs are resolved. 



3.2.1. Volume-limited unpaired field sample 

We construct a volume-limited sample of galax- 
ies with no detected companions (the "field" sam- 




Fig. 6. — Absolute magnitude distribution for 
galaxies in the full major pair sample (open) and 
the volume-limited major pair sample (filled). 



pie) for comparison with the volume-limited pair 
sample. The field sample includes only galaxies 
without companions that satisfy the pair selection 
criteria. The field galaxies in the volume-limited 
sample meet the same magnitude {Mr < —20.8) 
and redshift {z — 0.080 — 0.310) requirements as 
the major pair galaxies. The range of magni- 
tude differences excluded from the field galaxies 
{\AMr\ < 2) is slightly broader than the range ac- 
cepted for major pairs (|AMr| < 1.75) to exclude 
systems that could marginally be considered ma- 
jor pairs. We exclude galaxies within 70 kpc from 
the survey boundary or 70 kpc from the bound- 
aries around bright stars. The volume-limited field 
sample includes 2,234 galaxies. 

The distributions of Mr in the volume-limited 
major pair sample and volume-limited field sam- 
ple are similar. We use the Kolmogorov-Smirnov 
(K-S) test to determine that the two distribu- 
tions are drawn from the sample parent sample 
{Pks — 0.37). The volume-limited major pair 
sample and the unpaired field sample also have 
similar distributions of redshifts {Pks — 0.29). 
Figure [7] shows the redshift distribution of field 
galaxies. 

Some of the field galaxies may have a major 
companion fainter than the survey limit. We esti- 
mate the fraction of field galaxies likely to have 
an undetected major companion. The fraction 
of galaxies in pairs in the volume-limited sample 
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4.1. Identification of AGN 









1 




b 1 






J 





0.10 0.15 0.20 0.25 0.30 0.35 



Fig. 7. — Redshift distribution for the field galax- 
ies (open) and the volume-limited field galaxies 
(filled). 



and with apparent magnitude m R < 18.55 is 11% 
(102/902). The m R < 18.55 limit allows identifi- 
cation of pair galaxies with companions fainter by 
up to 1.75 mag. There are 1633 out of 2220 field 
galaxies (74%) in the range 18.55 < m,R < 20.3 
that could have a major companion fainter than 
our survey limit. An additional factor of 0.5 ac- 
counts for a pair galaxy having a 50% probability 
of being brighter than its companion. Hence the 
fraction of field galaxies likely to have an unde- 
tected major companion is ~ 4%. There are also 
some field galaxies in later stages of merging such 
that a companion is unresolved. Because the field 
sample contains some pairs and merger remnants, 
any differences in galaxy properties between the 
pairs and field samples attributable to interactions 
are lower limits. 

4. Spectroscopic Identifiers 

In this section we describe galaxy classification 
metrics to identify candidates for tidally triggered 
star formation. We identify the galaxies with 
strong AGN to remove them from our sample of 
star- forming galaxies ( §4.1|) . We describe the use 
of the spectroscopic indicator £>„4000 as a galaxy 
classification tool ( §4.2[) . 



We exclude strong AGN from our analysis be- 
cause measurements of emission lines, such as Ha, 
would not accurately represent star formation ac- 
tivity in galaxies with active nuclei. We distin- 
guish narrow-line AGN from star forming galaxies 
with measurements of the emission lines Ha, H/3, 
[N II] , and [O I I I] usin g the classification metric 
of iKewlev etai] (120061) . b ased on the "BPT" di- 
agrams of lBaldwin et all (|l98lh . The line fluxes 
Ha and Hf3 are stellar absorption corrected. 

The requirements we place on the spectra for 
classification are very mild: EW(Ha) > 2 A and 
EW(H/3) > 2 A. We impose no requirement on [N 
II] or [O III] because galaxies without measurable 
[N II] or [O III] could only have a very weak AGN. 
Furthermore, galaxies with EW(Ha) < 2 A have 
minimal star formation. Figure [8] shows the clas- 
sification diagram for the pair galaxies. 

We include composite galaxies in our pair and 
field samples, but exclude the galaxies classified 
as AGN. Although these criteria may admit weak 
AGNs, most of the contribution to the Ha emis- 
sion should be from star- forming regions. There 
are 119 narrow-line AGNs in the z = 0.080-0.376 
sample, 19 of which are in major pairs f ^3.2[) . 

Our ability to detect narrow-line AGN within 
the host galaxy declines with redshift because 
the covering fraction of the spectroscopic fiber in- 
creases with redshift. At redshift z = 0.080, the 
l'/5 fiber diameter corresponds to 2.2 kpc, and at 
z = 0.310, the l'/5 fiber covers 6.8 kpc. We also 
admit fewer AGN at higher redshift into our spec- 
troscopic sample because objects with stellar light 
profiles are excluded from our spectroscopic cata- 
log. 

Despite these caveats to our ability to detect 
narrow-line AGN at higher redshift, the mean red- 
shift of the narrow-line AGN in the full spectro- 
scopic sample (both pair and unpaired galaxies) 
is higher than that of the objects classified as 
starforming or composite: mean z — 0.27 and 
z = 0.23, respectively. The redshift distributions 
differ significantly (P KS = 1.6 x 10~ 5 ). We at- 
tribute this trend to the evolutiona ry decline in 
AGN activity at low redshift (e.g. iBrand et al. 
20051 : ICavalierell2000h . 

We identify galaxies with broad emission lines 
by visually inspecting all spectra that are potential 
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Slorforming ± 
Composlte A 
O Norrow-line AGN<> 




-1.0 -0.5 0.0 0.5 
Log ([N ll]/H«) 



Fig. 8. — Classification of pair galaxies as star- 
forming, composite, or narrow-line AGN using the 
Kewlev et al.l ( 20061 ) criteria. The solid black line 



shows the boundary between composite galaxies 
and AGN, and the dotted line shows the bound- 
ary between star forming and composite galaxies. 
We include composite galaxies in our sample for 
analysis of star formation activity. 



broadened line objects. The criteria for broad-line 
AGN candidates is EW(Hq) > 3 A and the line 
flux in a band widened from ±8 A to ±12 A in- 
creases the measured Ha flux by more than 10%, 
or EW([N II]) > 3 A and the line flux in a 
band widened from ±8 A to ±12 A increases 
the measured [N II] flux by more than 10%. We 
identify 26 galaxies with broadened lines in the 
z = 0.080 — 0.376 sample. Five of these are in 
major pairs. 

We exclude galaxies classified as AGN from our 
analysis of star formation activity; we retain the 
non-AGN companion. Table [3] reports the AGN 
fraction in each pair sample and in the matching 
unpaired field sample. The AGN fraction in pairs 
exceeds that in the matching field sample for each 
of the subsamples. The increase in AGN fraction 
for both the full and volume-limited major pair 
samples is a factor of ~ 2, a significance of ~ 2a. 
The higher AGN fractions in the volume-limited 
samples (pair and field) compared to the full sam- 
ples is consistent with the exclusion of lower lumi- 
nosity galaxies from the volume-limited sample, 
which are less likely to contain an AGN. 



4.2. Galaxy classification by D„4000 

It is useful to segregate galaxies into two cate- 
gories: early-type galaxies, which are generally 
gas-poor and have little or no active star for- 
mation, and late-type, which contain more gas 
and have young stellar populations. The gas- 
rich systems are potentially susceptible to tidally 
trig gered star formation in major interactions 



Mihos fc Hernquistlll996l;lTissera et al.ll2002l 

20071) . 



Cox et all 120061: iDi Matteo et al.l 120071 ). Inter- 
actions between gas-poor galaxies ("dry merg- 
ers") produce little or no star formation activ- 
ity, although they contribute su bstantially to the 

TVan et alj|2005t 



build-up of massive galaxies (e.g. 
IVan Dokkuml 120051 : 



Cattaneo et al 



2008). 



Our 

analysis of star formation activity in pair galaxies 
focuses on the late-type systems. 

There are a number of classification schemes 
to separate the early and late-type galaxy pop- 
ulations. Photometric discriminants include 
color, concentr a tion, a nd absolute magnitude (e.g . 



IStrateva et al.l l200lfc iKauffmann et al.l l2003bl ). 



and spectroscop ic methods include P„,4 000 and 
H<5 absorption ( Kauffmann et all 12003a ). The 
Z)„4000 indicator discriminates by stellar popu- 
lation age. At wavelengths bluer than 4000 A, 
metal lines in low mass stars absorb the light and 
cause a "break" in the spectrum. As the stellar 
population ages and the massive, hot stars die 
off, D ra 4000 incr e ases m onotonically with time. 



Kauffmann et al.l (|2003al ) use stellar population 
models to show that galaxies with Z)„4000 < 1.5 
have young stellar populations (< 1 Gyr). Metal- 
licity has a strong effect on the the value of £>„4000 
only after 1 Gyr past a burst of sta r formation (see 



Fig. 2 in IKauffmann et al.ll2003al ). Measurement 



of D„4000 is insensitive to galaxy reddening. 
(l2008h 



Vergani et al 



use D„4000 to separate 
spectroscopic early-type galaxies from late-type 
galaxies at a dividing line of D„4000 = 1.5 in their 
analysis of galaxy stellar mass assembly in th e VI- 
MOS VLT Deep Survey. iMignoli et all(|2005l) also 
use ZM000 = 1.6 as a dividing line (different def- 
inition of D4000 from lBruzualHl983T : ratio of flux 
in bands 4050-4250 A to 3750-3950 A) along with 
other spectral measurements to classify galaxies in 
the K20 survey, a near-IR selected redshift survey. 
We follow this approach. 

The distribution of £>„4000 is bimodal, sepa- 
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Table 3 
agn fraction 



Sample 


Pair a 


Field b 


Significance 


Pair sample size 


Full 


8 ±2% 


4 ±2% 


2a 


280 


Volume limited 


12 ± 3% 


5 ±2% 


2a 


128 



a Fraction of AGNs compared to fraction of galaxies that meet our 
criteria for classification f ^4.ip . We identify only very bright AGNs with 
our spectra and our classification criteria. 

b Proportional representation of field galaxies matching the Mr and 
redshift distribution of the pair sample. 



rating galaxy populations dominated by old stars 
from systems with recent star formation. Figure^ 
shows the distribution of £>„4000 for the 6,644 
galaxies (both pair and unpaired) at z = 0.080 — 
0.376 with a robust D„4000 m easurement. The 
peaks a re at ~ 1.15 and ~ 1.75. IKauffmann et al 



( 2003al ) similarly observe a bimodal distribution in 
D„4000, with peaks at 1.30 and 1.85 in their sam- 
ple of - 100, 000 galaxies in the SDSS. We choose 
the minimum between the bimodal distribution as 
our dividing line, D„4000 = 1.44 (Figure [9]). We 
refer to galaxies with Z?„4000 < 1.44 as "low" 
D„4000 galaxies and those with L>„4000 > 1-44 
as "high" Z?„4000 galaxies. Our analysis of star 
formation activity includes the low Z?„4000 galaxy 
in a mixed pair, but does not require both galaxies 
to have D n 4000 < 1.44. 

The emission line fraction is a strong function 
of D„4000. Figure [10] shows the steep decline 
in the fraction of emission line galaxies between 
D„4000 = 1.3 — 1.5. The emission line galaxies 
have EW(Ha) > 3 or EW([0 II] )> 3. Thus segre- 
gating by Z?„4000 is reasonable and corresponds 
well to segregating by the presence of emission 
lines. 

We also compare the use of -D„4000 for galaxy 
classification with color, another widely used in- 
dicator of galaxy type. Figure [11] shows the bi- 
modal distribution of rest-frame SDSS (g — r) 
color versus D„4000. We compute the rest-frame 
SDSS (g — r) color using the k+e correction de- 
termined by lA nnis (2001) from the Pegase code 
(|Le Borgne fc Rocca-Volmeragell2002l ). The sep- 
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Fig. 9. — Distribution of £>„4000 for all galaxies 
in the spectroscopic sample. The bimodal distri- 
bution has a local minimum at £>„4000 = 1.44 
(dashed line). 



aration by D„4000 is more sharply defined than 
that of rest-frame color. The rest-frame galaxy 
colors are affected by reddening, and depend on 
the noisy k+e corrections. Z?„4000 provides a well- 
defined method to segregate galaxies and has the 
advantage that no redshift dependent corrections 
are required. 

5. Local Density Effects 

Because pair galaxi es are more comm on in 



higher density regions ([Barton et all l2007h . and 



because cluster and group galaxies are more likely 
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Fig. 10. — Fraction of emission-line galaxies as a 
function of D„4000 for all galaxies in the spec- 
troscopic sample, excluding AGN. Emission line 
galaxies have EW(Ha) > 3 or EW([0 II] )> 3. 
There is a steep decline in the fraction between 
-D„4000 = 1.3 — 1.5. Error bars from bootstrap 
resampling indicate the 95% confidence intervals. 



to be red with less active star formation than 
isolated galaxies (e . g^ 



Hubble fc Humason 1931 



Cooper et al.l 120071 : iGerke et al.l 120071 ). compar 



isons between pair and field populations without 
attention to density can strongly b ias the inter- 
pretat ion of studies of close pairs. iBarton et al 



(|2007l ) emphasize that including high density re- 
gions suppresses the detection of triggered star 
formation. 

To compare the environment of pair and field 
galaxies, we compute the number of companions 
(N c ) within a co-moving sphere of radius 985 kpc 
centered on the galaxy (radius equivalent to 700 
/ij^gkpc). The 985 kpc radius is withi n the typical 



virial radius of clusters ( 
2003), and maximizes the survey area included in 



1 /t 10 1 Mpc. lRines et al 



our analysis. This density measureme nt is consis- 



tent with that of IBarton et al.l (|2007l ). We mea- 



sure N c within the volume-limited sample, requir- 
ing that the galaxy reside within z = 0.080 — 0.310 
and have magnitude Mr < —20.8. We count all 
neighboring galaxies, not just major companions. 
We exclude regions within 985 kpc of the survey 
edge. 



Barton et alJ ((2007) predict that most observed 



pair galaxies reside in higher mass cluster or 



group-sized halos, whereas field galaxies are usu- 
ally in isolated low mass halos. We test this pre- 
diction with our complete volume-limited data set, 
using N c as a proxy for halo membership. Note 
that our N c counts neighbors; the halo count N 
in Barton et al. counts total halo occupation, i.e. 
N = N c + 1. 

Figure [12] shows the fraction of galaxies with 
various N c as a function of projected distance to 
the nearest neighbor. Pair galaxies fall in the 
range AD < 70 kpc. As Barton et al.l show, pairs 
occur more frequently in locally dense regions. 
We find that 32% of pair galaxies have N c > 8. 
The overall fraction of galaxies in dense regions 
is smaller: 19% of all galaxies lie in regions with 
N c >8. 

Our observations are in excellent agr eement 
with the predictions of Barton et al. ( 20071 ). They 
find that 39% of pair galaxies are in a host halo 
with total number of galaxies N > 9; 19% of galax- 
ies in the sample as a whole have N > 9 (N c > 8). 
Selecting a sample of field galaxies for a fair com- 
parison of specific star formation rates thus neces- 
sitates consideration of local density effects. 

Selecting star forming galaxies with young stel- 
lar populations (low £>„4000) already reduces the 
effects of the environment because young, blue 
galaxies are relatively more abundant in low den- 
sity regions. In the volume-limited major pair 
sample, 63% of the low £>„4000 non-AGN galaxies 
reside in regions with N c < 4. 

Restricting the analysis of star formation ac- 
tivity solely by -D„4000 does not, however, elim- 
inate differences in star formation activity with 
local density. The distribution of local densities 
differs between low D„4000 pair and field galaxies 
(Pks = 1-6 x 10~ 5 ). Thus we restrict any direct 
comparisons of specific star formation activity in 
pairs and field galaxies to low density regions to 
ensure that density effects do not dominate our 
measurements of star formation activity. A mod- 
erate limit of N c < 4 maximizes our sample size. 

6. Measurements of tidally triggered star 
formation 

Here we quantify the frequency, strength, and 
timescale of the triggered star formation for galax- 
ies at intermediate redshift. We examine trends 
in star formation activity across redshift and lu- 
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mmosity for the galaxies in the volume-limited 
sample. The completeness and quality of spec- 
trophotometry enable us to carry out the most 
detailed spectroscopic analysis of pair galaxies to 
date in the redshift range z — 0.080 — 0.376. 
We examine the Ha specific star formation rate 
(SSFRj/ Q , §2.3. lj) . the spectroscopic parameter 
D„4000 ( §4.2j) . and a set of stellar population 
models ( f2XT]l . 

The volume-limited pair sample is well suited 
to comparison with the predictions of numerical 
simulations. We c ompare our obse r vation s with 



the predictions of Di Matteo et al. ( 20081 ) . who 



measure the intensity, frequency, and duration of 
merger-driven star formation in their large suites 
of numerical simulations of major interactions. 

6.1. Star formation indicators and pro- 
jected separation 

We examine a set of star formation indica- 
tors as a function of projected separation to look 
for the signature of triggered star formation at 
redshifts 0.080 — 0.376. The star formation - 
AD anti -correlation frequently observed at low 
redshift jBarton et all l2000t iLambas et all [2 003; 



Nikohc et al.ll2004UWqods et al.l2006HGeller et al 



20061: IWoods fc Gellerll2007l : lLi et al. 2008a; Ellis on et al.l 1.1, iLin et all (|2007l ) 



2008) results from an increase in central star for- 
mation activity triggered by a close pass from a 
neighboring galaxy. As the pair galaxies move 
apart and the burst ages, the star formation ac- 
tivity decreases. 

We observe a strong anti-correlation between 
SSFR^a and AD in the sample of all low D„4000 
major pair galaxies at z = 0.080 — 0.376 (at all lo- 
cal densities). The Spearman-rank test computes 
a probability of no correlation of Psb. = 6.0 x 10~ 4 
for the 134 galaxies in the sample (Figure [T3l top 
left panel.) The range of absolute luminosities and 
local densities in this sample corresponds to that 
included in typica l samples at low redshift (e.g. 
Barton et allboOOl ). 

Figure [13] shows the mean SSFR# Q as a func- 
tion of AD for low D„4000 galaxies in low density 
regions in the volume-limited major pairs sample 
(top right panel). We measure a correlation be- 
tween SSFR# Q and AD for the 70 galaxies in this 
sample, where Psb — 1.3 x 10~ 2 . The decreased 
significance compared to the sample of all major 



pair galaxies results from the reduction in sam- 
ple size. Excluding the lowest luminosity galaxies 
from the volume-limited sample may also reduce 
the signal of the interaction because low luminos- 
ity galaxies are more strongly affected than the 
more luminous galaxies. 

We measure a weak correlation between AD 
and D„4000 for the low D„4000 galaxies in the 
volume- limited pair sample at the level of P$ b. = 
5.1 x 10~ 2 ; bursts with small AD have small 
D„4000 (Figure[13l bottom left panel). This trend 
is consistent with the expectation that close pairs 
have had the most recent starbursts, and hence 
harbor the youngest stellar populations. 

We apply another measure of the recent star 
formation based on the stellar population mod- 
els. The model-determined fraction of total lu- 
minosity at 5500 A from the youngest discrete 
starburst population age (5 Myr) correlates with 
AD for the volume-limited major pair sample: 
P SR = 1.4 x 10~ 2 . Pairs at the smallest AD 
have the highest fractional contribution from a 
5 Myr stellar population, compared to galaxies 
with larger AD (Figure [T3l bottom right panel). 

These results are consistent with other re- 
cent observations. At redshifts of 0.1 < z < 
likewise observe an anti- 
correlation between the median infrared lumi- 
nosity of merging galaxies and pair separation in 
their sample of ~ 100 systems with data from the 
DEEP2 Galaxy Redshif t Surv ey and HST/ACS 
imaging. Ide Ravel et~ai1 (|2009h similarly find that 
the galaxy pairs at the smallest separations have 
the greatest median EW([0 II]) in their sample 
of 251 pairs in the VIMOS VLT Deep Survey 
( Le Fevre et al. 20051) . which has a mean of red- 
shift ofz = 0.76. 

We next explore properties of the galaxies that 
drive the correlation between star formation in- 
dicators and AD — the galaxies with significant 
star formation rates at small separation. The 
galaxies with AD < 25 kpc and SSFR# Q > 
4 M Q yr _1 /10 10 M Q have luminosities similar to 
the rest of the low D„4000 galaxies in the volume- 
limited sample (Mb ~ —21.4). These galaxies are 



bright; Mt 



-22.1 for all galaxies in our spec- 



troscopic sample in this redshift range (Diaferio 
et al. 2010, in prep.). The correlation between 
SSFR# Q and AD for luminous galaxies at red- 
shifts z = 0.1 — 0.3 is an interesting observation 
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for comparison with lower redshift results. 

In the low redshift CfA2 Redshift survey, the 
lower mass galaxies (measured by rotation veloc- 
ities) exhibit a anti-correlation between E W(Hoi) 
and A D; more massive galaxies do not ((Barton 
(|l999h ). They conclude that the least massive 
galaxies exhibit the strongest bursts of star forma- 
tion and are responsible for drivi ng their observed 
EW(H a)-A £> correlation . Both IWoods fc GeUer 
(|2007t) and lEllison et all (|2008h find at low red- 
shift that low mass (luminosity) galaxies exhibit 
relatively more powerful triggered specific star for- 
mation than high mass galaxies. 

Finding luminous galaxies with strong evi- 
dence of triggered star formation at interme- 
diate redshift is consistent with "downsizing" 
(|Cowie et al. l ll996h . which suggests that higher 
mass (luminosity) galaxies formed their stars 
earlier, and that the lower mass galaxies have 
more efficient star formatio n at later time s than 
the higher mass galaxies (iGuzman et al 



Brinchmann k Ellis! l2000t iKodama et al 



1997 



2004: 



Bell et all 120051 : Uuneau et all 120051: iNoeske et al 
2007). 



6.2. Frequency and strength of triggered 
star formation 

Gravitational interactions clearly trigger star 
formation in some cases (see references in fjl}. 
Some galaxies fail to have enhanced star for- 
mation because the pairs have not yet reached 
perigalacticon, the pairs may be mere interlop- 
ers along the line-of-sight, or they may have in- 
ternal structure less cond ucive to gaseous inflows 
(|Mihos fc Hernauistl[l996l) . 

We place limits on the frequency of triggered 
star formation by comparing the SSFR# Q of the 
low D„4000 pair galaxies with those of the low 
D„4000 unpaired field galaxies in the volume- 
limited samples. Figure [14] shows the normal- 
ized distributions of SSFR# Q for the low D„4000 
galaxies in major pairs and the unpaired field 
galaxies. The dotted line indicates the median 
SSFR# Q for the field galaxies. Table 2] lists the 
quartiles of the distribution and the fraction of 
pair galaxies with SSFR# Q more than twice the 
field median. The results for the very close pairs, 
AD < 25 kpc, are listed separately. 

Table 3] shows that the low D„4000 major pair 



galaxies have an excess of high SSFR# Q compared 
to the low D„4000 field galaxies. We define the 
fraction of pair galaxies experiencing enhanced 
specific star formation as SSFR# Q = 2 x F50, 
where F50 is the median SSFR# Q for the field 
galaxies. According to this definition, 32 ± 7% 
of major pair galaxies experience enhanced spe- 
cific star formation. A greater fraction of close 
pairs with AD < 25 kpc, exhibit enhanced spe- 
cific star formation rates: 42 ± 13%. Observing a 
greater fraction of close pairs with enhanced star 
formation is consistent with triggering at the clos- 
est approach and the subsequent decline as the 
burst ages and the pair moves apart. We find only 
one galaxy with SSFR# Q > 5 x F 50 - 

A number of generic selection issues affect the 
measured frequency and strength of triggered star 
formation. The first issue is that a small fraction 
(< 4%) of the field galaxies is likely to be in a 
pair with an undetected companion fi j3.2.1[) . Sec- 
ond, excluding galaxies with strong AGN preferen- 
tially excludes pairs with active star formation be- 
cause triggered AGN activity and star formation 
activ ity often occur as p art of the same process 



(e.g. Hopkins et al. 20081) . Third, we are unable 



to resolve some pairs or mergers at separations 
AD < 15 kpc, where the stro ngest star formatio n 
enhancement is expected (e.g. lBarton et ai] l2000). 
Our measurement of the frequency and strength of 
triggered star formation should therefore be inter- 
preted as lower limits. 

We compare measurement of the frequency 
of triggered star formation from our volume- 
limited major pa i rs sam ple with the predictions of 
Di Matteo et al.l ( 2008 '). who use numerical sim- 



ulations to study the frequency, intensity, and 
duration of triggered star formation activity. Be- 
tween 25 and 50% of their "fly-bys" have star 
formation rates twice the isolated case (see their 
Table D.l.) We compare with their fly-bys and 
not their merger scenario because the galaxies in 
our sample are distinct systems, which have not 
yet had a final merger. Our results are consistent 
with the predictions of Di Matteo et al. 

One difference between our observations and 
the Di Matteo et al. simulations is that their max- 
imum star formation rate refers to the lifetime of 
the galaxy; ours represents a snapshot in time. 
Because the most intense bursts of star formation 
occur over relatively short time scales in the course 
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Table 4 

Distribution of specific star formation rates in pairs and unpaired field galaxies 



Sample 3, 


Q 25 b 


Q50 


Q75 


> 


2F 50 C 


Sample size 


Pair 


1.55 


2.61 


4.26 


32 


±7% 


72 


Close pair d 


1.87 


3.11 


4.43 


42 


±13% 


26 


Field 


1.18 


1.80 


2.86 


14 


±1% 


777 



a Pair and field samples derive from the volume-limited sam- 
ple of low £>„4000 galaxies, where Mr < —20.8. Pair galaxies 
are in major pairs, |AMr| < 1.75. 

b Quartile of the distribution of SSFR# Q in units of 
M Q yr- 1 /10 w M G . 

c Fraction of pair galaxies with SSFRff Q > 2F50, where F50 
is the median of the field galaxies in the same redshift bin. 

d Close pair: AD < 25 kpc. 




Fig. 11. — D„4000 versus rest- frame (g — r) color. 
Two galaxy populations are visible: systems with 
blue colors (lower left), and systems with red 
colors (upper right). The dashed line indicates 
£> n 4000 = 1.44 (see Figure E). 



of the merger (see Figure 4 in Di Matteo et al.), we 
are unlikely to observe pair galaxies at maximum 
intensity. We therefore expect to measure a lower 
frequency than that predicted by Di Matteo et al. 
This prediction is consistent with our finding only 
one galaxy with SSFR^ Q > 5 x F 50 ; Di Matteo et 
al. find that 15% of their major mergers have star 
formation enhanced by this large factor. 

The frequency and strength of triggered star 
formation that we measure are consistent with the 



observations of iJogee et al.l (|20081 ). who find that 
the average star formation rate of strongly dis- 
turbed galaxies exhibits only a modest increase 
over the morphologically undisturbed galaxies in 
their sample of ~ 4, 500 galaxies at 0.24 < z < 
0.80. Our results are also in line with the low 



redshift observations of lEllison et al.l (|20081 ). who 
measure star formation rate enhanceme nt in SDSS 
galaxies at z < 0.16. Li et al. I (12008 a) find that 
the SDSS systems with the highest star formation 
rates are likely to have a companion, but not all 
systems with close companions have high star for- 
mation rates. 

A consistent explanation for the range of re- 
sults drawn from simulations and observations is 
that strongly enhanced star formation is rare and 
short lived. Simulations can track the maximum 
enhancement across the lifetime of the merger; dif- 
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ferent types of pair selection probes systems in dif- 
ferent stages of the interaction. Density effects 
may also reduce th e impact of observat ions of star 
formation in pairs (jBarton et al.ll2007l ). 



6.3. Duration of triggered star formation 

We use the stellar population models described 
in ij2.3.1l to compare the stellar composition of pair 
and unpaired field galaxies. The stellar population 
models fit the spectra with a discrete set of star- 
bursts of age 0.005, 0.025, 0.1, 0.3, 0.6, 0.9, 1.4, 
2.5, 5 and 10 Gyr. From the models we extract 
the contribution of each starburst population to 
the flux at 5500 A. 

Figure [15] shows the ratio (pair/field) of the 
mean fraction of flux at 5500 A attributed to each 
discrete starburst as a function of starburst age for 
the starburst populations included in the model. 
The pair galaxies clearly contain a larger frac- 
tion of young stellar populations up to burst ages 
~ 300 — 400 Myr. The ratio of the mean fraction of 
flux from each starburst population dips slightly 
below one for burst ages > 500 Myr, because the 
younger stellar populations contain a greater frac- 
tion of the flux in the pair galaxies. 

Our measuremen t agrees well with the results of 
Barton et al.l (12000), who app l y spectral synthesis 



models (Lcithcrcr et al. 1999; Bruzual & Chariot 



1999) to their data fr om the CfA2 Redshift Sur- 
vey (|Falco et al. 1999). They determine that the 
Ha emission and the galaxy colors are best de- 
scribed by a burst with a continuous duration 
> 100 Myr on top of the p re-existing stellar popu- 



lation. The simulations of Di Matteo et al. (2008) 



similarly suggest a merger-driven starburst dura- 
tion of up to a few hundred Myr, consistent with 
our measured duration. 

7. Summary and conclusion 

We examine spectroscopic properties of pair 
and field galaxies to quantify effects of gravita- 
tional interactions at intermediate redshift. Our 
sample derives from th e Smithsonian Hectospe c 
Lensing Survey fSHELS: lc"eTler et al.l[2005l . l2010l) . 
SHELS includes 9,825 galaxies and is 97.7% spec- 
troscopically complete to R = 20.3 over an area 
of 4 deg 2 . We select for galaxies in the redshift 
range z = 0.080 — 0.376. This survey represents 
the most complete spectroscopic sample in its red- 



shift range. 

We focus on the systems that have the poten- 
tial to exhibit bursts of star formation as a re- 
sult of the interaction. Substantial evidence shows 
that major pairs are more strongly affected by 
the interaction. We identify a full set of ma- 
jor (|AMr| < 1.75) pairs including 622 galaxies 
in the redshift range z = 0.080 — 0.376, and a 
volume-limited subset of major pairs in the red- 
shift range z = 0.080 — 0.310, including 327 galax- 
ies to M R = -20.8. 

Within our major pair sample, we further nar- 
row our selection using the spectroscopic index 
£>„4000 = 1.44 as the divide between systems with 
older stellar populations, and systems with young 
stellar populations that likely contain gas. We fur- 
ther restrict our sample to systems with low sur- 
rounding density, which we measure with a count 
of neighbors within a volume of comoving radius 
985 kpc. 

The spectroscopic diagnostics of Ha specific 
star formation rate (SSFR# Q ), D„4000, and a set 
of stellar population models enable the investiga- 
tion of the strength, frequency, and timescale of 
triggered star formation. We show: 

1. The spectroscopic indicator Z) n 4000 pro- 
vides a useful classification metric and corre- 
sponds closely with identification of emission 
line galaxies. 

2. The star formation indicators SSFR# Q , 
Z?„4000, and presence of young stellar popu- 
lations exhibit an anti-correlation with AD, 
demonstrating that bursts of star formation 
are associated with close proximity to a ma- 
jor companion. 

3. 32 ± 7% of major pair galaxies in the volume- 
limited sample experience enhanced specific 
star formation activity at twice the median 
of the unpaired field galaxies. For very close 
pairs (AD < 25 kpc), the fraction is 42 ± 
13%. This trend is consistent with the tidal 
triggering picture. 

4. We use stellar population models to show the 
burst of star formation following an interac- 
tion has a duration of ~ 300 — 400 Myr. Pair 
galaxies show an increase over the field in the 
light fraction from young stellar populations 
for burst ages up to ~ 300 — 400 Myr. 
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The most effective way to increase our ability to 
measure differences between pair and field galax- 
ies as a function of redshift, or to determine the 
AGN fraction in pairs would be to observe a larger 
population of very close pairs (AD < 15 kpc) at 
redshift z ~ 0.3 using small aperture spectroscopy. 
It is important to have high resolution photomet- 
ric data in combination with good seeing to dis- 
tinguish close pairs. 
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Fig. 12. — Fraction of galaxies in the volume- 
limited sample with a total iV c neighboring galax- 
ies within 985 kpc co-moving (equivalent to 
700 h^ kpc) and AV < 1000 km s^ 1 for the 
volume-limited sample as a function of distance 
to the nearest neighbor. Pair galaxies have AD < 
70 kpc. 
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Fig. 15. — Ratio of mean light fraction from star- 
burst populations as a function of starburst pop- 
ulation age for pair and unpaired field galaxies 
with low £>„4000 in the volume-limited sample. 
The pair galaxies show a greater light fraction 
from young stellar populations up to burst ages 
<~ 300 — 400 Myr. The dashed line indicates a 
ratio of unity. Error bars are from bootstrap re- 
sampling: boxes indicate the inter-quartile range 
and the outer lines show the 95% confidence inter- 
val. 
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Fig. 13. — Mean star formation indicators versus projected separation. We show the SSFR# Q versus pro- 
jected separation for low D„4000 major pair galaxies in the full sample at all local densities (top left). Star 
formation indicators for the low D„4000 galaxies in regions of low local density in the volume-limited major 
pair sample: SSFR^ Q (top right), £>„4000 (bottom left), and fraction of light at 5500 A from a 5 Myr star- 
burst population (bottom right). Error bars are from bootstrap re-sampling: boxes indicate the inter-quartile 
range and the outer lines show 95% confidence intervals. For each star formation indicator shown above, the 
galaxies exhibit a correlation between the star formation indicator and AD. 
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Fig. 14.— Normalized distribution of SSFR Hq for the low L>„4000 pair (N = 72) and field galaxies (N = 
777) in the volume-limited sample. The pair galaxies' distributions exhibit a clear excess of high SSFR# Q 
compared to that of the field galaxies. The dotted line on both plots indicates the median SSFRhq for 
the field galaxies. The dashed line on the pair galaxies' plot shows the contribution from galaxies with 
AD < 25 kpc. 
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